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The branch of analytic number theory which is concerned with the number of integer points in large 
domains, in the Euclidean plane as well as in space of dimensions > 3, has a long and very prolific 
history which reaches back to classic works of E. Landau, J.G. Van der Corput, G. Vorono'i, G.H. Hardy, 
to mention just a few. Enlightening accounts of the developments of this theory can be found in the 
monographs of F. Pricker PH, E. Kratzel [EH, [EH, and M. Huxley g2|. Although the major problems in 
this field are quite old (and yet unsolved), there have been a lot of amazing new developments in recent 
times. This survey article attempts to give an overview of the state-of-art in this field, with an emphasis 
on results established during the last few years. Furthermore, we are aiming to point out the intrinsic 
relationships between lattice point quantities and certain arithmetic functions which arose originally from 
number theory originally without any geometric motivation. 



1. Lattice points in a large circle. The problem of Gaufi. 

Starting from the arithmetic function viewpoint, let 

r(n) := #{(m, k) £ 1? : m 2 + k 2 = n } 

denote the number of ways to represent the integer n > as a sum of two squares of integers. While 
the maximal size of r(n) is estimated, e.g., in E. Kratzel |79j. and an explicit formula may be found in 
E. Hlawka, J. Schoifiengeier, and R. Taschner [3^1, p. 60, its average order is described by the Dirichlet 
summatory function 

A{x) = r ( n ) > 

0<n<x 

where x is a large real variable. In geometric terms, A(x) is the number of lattice points (m, k) £ Z 2 
in a compact, origin-centered circular disc C(y/x) of radius \/x. Obviously, A(x) equals asymptotically 
the area of C(y/x), i.e., irx. Already C.F. Gaufi (1777 - 1855) made the following more precise but still 
extremely simple observation: If we associate to each lattice point (m, k) £ Z 2 the square 

Q(m, k) = {(u, v) £ R 2 : \u-m\<\, \v-k\<\}, 

then obviously 

C 



(m,fc)eC(VH) 

Evaluating the areas, we see that, for every x > i, 



VS-^)c |J Q(m^)cc(vS+|] 



hence the lattice point discrepancy of C(^/x) (traditionally often mis-translated as "lattice rest") 

P(x) :— A(x) — irx (1.1) 

satisfies 

\P(x)\ <irV2x~+^ (1.2) 

which is usually simply written as P(x) = O(yfx) or P{x) <C ^fx. 

Before reporting on stronger results, we shall try to shed some light on the mathematical background 
of all more sophisticated approaches to the problem. The naive idea is to apply Poisson's formula to 
evaluate the sum 

A{x) = ^ c x(m, k) , 

(m,fc)£Z 2 

where c x denotes the indicator function of C(y/x). To get rid of the discontinuity at the boundary of 
C(y r x), which would create subtle questions of convergence, it is convenient to consider the integrated 
form 

X X / \ 

I A(t)dt= f I J2 1 ) dt = J2 (x-m 2 -k 2 ) = 







m 2 +/c 2 <i / rn 2 -\-k 2 <x 



= H J J (x-u 2 -v 2 )e 2 ^ mu+kv U(u,v). 

(m,k)£Z 2 u 2 +v 2< x 

Since the double integral can be evaluated explicitly by means of the Bessel function J2, this gives 

/ A(t) dt = - x 2 + - Y ^ J 2 (2n^) . (1.3) 
J 2 7r ^— ' n 

n=1 

With some care, one can verify that this identity may be differentiated term by term (cf., e.g., E. Kratzel 
|83| . p. 126, Th. 3.12), yielding (for x not an integer) 

00 

P(x) = yj^V(n)n- 1/2 Ji(27iV^) , (1.4) 



71=1 



which is called Hardy's identity. Since this latter series is only conditionally convergent, it is convenient 
for applications to have a sharp truncated version at hand. A result of this kind can be found as Lemma 
1 in A. Ivic |55] : For x > 1, x ^ Z, and x < M < x A , A > 1 some fixed constant, 

P(x) = y/x r{n)n- 1/2 J 1 {2TTy/nlE) + 

l<n<M 

+0 (min (V^M- 1 / 2 + x^+'M-^Wxir 1 + z^M" 1 / 4 , x e )) , (1.5) 

where || • || denotes the distance from the nearest integer, and e > is arbitrary. Moreover, using the 
well-known Bessel function asymptotics 

Ji(27rV^F) = -(nx)~ 1/4 cos(27iV«x- §tt) + O ({n X y 3/i 
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one arrives at a representation 

P(x) = — x 1 / 4 ^ ~Yfi cos(2ny/nx — |?r) + remainder terms . (1.6) 

n l<n<M 71 

In a sense, all of the deep investigations on the order of P{x) which have been done during the past 
century, are based on results like (1-6), though often somewhat shrouded by technicalities. 

1.1. O-estimates in the circle problem. For instance, to derive upper bounds for this lattice point 
discrepancy, one is led to consider exponential sums 

^(m 2 + k 2 r 3/4 e 2 ™^^+^ . 

m,k 

The twentieth century saw the development of increasingly complicated methods to deal with such 
sums. A historical survey of the results obtained can be found in E. Kratzel They started with 

W. Sierpihski's |156| bound 

P(x) = O (V /3 ) 

and in the 1980's culminated in the work of G.A. Kolesnik 74 , who had 

P(x) = O (V39/429+^ 

as his sharpest result. 

During the last one-and-a-half decades, M. Huxley, starting from the ideas due to E. Bombieri - H. 
Iwaniec jS] and ideas due to H. Iwaniec and C.J. Mozzochi devised a substantially new approach 
which he called the Discrete Hardy- Littlewood method. His strongest result was published in 2003 |43) 

P(x) =o(x 131 /416 (log:r) 18637/8320^ 

thereby improving upon his 1993 gT| bound 0{x 23 ^ 73+e ). Note that |f| = 0.3149038 . . ., §§ = 0.315068 . . ., 
while Kolesnik's ±§§ = 0.324009 .... 

1.2. Lower bounds for the lattice point discrepancy of a circle. As (1.6) suggests, a stronger 
estimate than P(x) — 0(x 1 ^ 4 ) cannot be true. In fact, G.H. Hardy [3U|, [3T] proved that 1 

P{x) = n ± (x 1/4 ), P(x) = n_{x 1/4 {\ogx) 1/4 ), (1.8) 

adding the unsubstantiated claim that P(x) = ^(^^(logx) 1 / 4 ). The second part of (1.8) was not 
improved until more than six decades later by J.L. Hafner |28j . who obtained 2 

P(x) = ft- (x l / 4 {\ogx) 1/4 {\og 2 x) (log2)/4 exp(-c(log 3 x) 1/2 )) , c> . (1.9) 

J Let us recall briefly the different C-symbols: For real functions F and G > 0, we mean by F(x) = Q^.(G(x)) that 
limsup(F(a;)/G(a;)) > 0, as x -> oo. Similarly, F(x) = f2_(G(x)) stands for -F(x) = Q + {G{x)). Further, F(x) = Q ± (G(x)) 
says that both of these assertions are true, and F(x) = Q(G(x)) means that at least one of them applies, i.e., F(x) ^ o(G(x)). 

2 We write throughout logj for the j-fold iterated natural logarithm. Thus log 2 = log log, and so on. 
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All estimates like (1.8) and (1.9) start from (a variant of) formula (1.6) and are based on the idea of 
finding an unbounded sequence of a;- values for which the terms cos(2iryJnx — §7r) "essentially all pull 
into the same direction". That is, the numbers y/nx should be close to an integer, which is achieved 
by Dirichlet's theorem on Diophantine approximation (see e.g., Lemma 9.1 of 60 ). In Hafner's work, 
this task is done only for those n for which r(n) is comparatively large; this idea slightly improves the 
estimate. 

As for r2 + -results, K. Corradi and I. Katai ^2] proved that 

P(x) =n + (a; 1 / 4 exp(c'(log 2 a;) 1 / 4 (log3a;)- 3 / 4 )) (c' > 0) , (1.10) 

thereby refining previous work of A.E. Ingham £HJ and of K.S. Gangadharan [251 - (Note that the factor of 
x 1 / 4 here is less than any power of logx. This inherent asymmetry arises from the fact that the deduction 
of (1.10) is based on sort of a quantitative version of Kronecker's approximation theorem which is much 
weaker than Dirichlet's.) 

The strongest f2-bound for P{x) available to date was established in 2003 by K. Soundararajan |157j . 
His ingenious new approach allows us to restrict the application of Dirichlet's theorem to a still smaller 
set of integers n. His method can be summarized in the following fairly general statement. 

Soundararajan's Lemma |157| . Let (f(n))'%L 1 and (X n )^Li be sequences of nonnegative real numbers, 
(X n )^—i non- decreasing, and 53 n =l f( n ) < 00 ■ Let L > 2 and Y be integer parameters, 7 S K fixed. Put 

00 

S(t) = f( n ) cos(27rA„t + 7) . 

n=l 

Suppose further that M is a finite set of positive integers, such that {X m : m S M. } C [jjAy,|Ay]. 
Then, for any real T > 2, there exists some t £ [\T, (6L)'^ V1 ' +1 T] such that 

\s(t)\ > - e / M - j— j y: /w-^e/w- 

m£M n: A„<2Ay n=l 

Furthermore, if '7 = wtt , ioeZ, then on the left hand side \S(t)\ can be replaced by (—1)™ S(t). 

When applied to (1.6), this lemma yields 

P(x) = Q[x 1 /\logx) 1 /\log 2 x)^ 2l/3 - 1 \log 3 x)- 5 / 8 ) . (1.11) 

To assess the refinement in the exponent of the log 2 x-factor, note that i log 2 = 0.1732 . . ., while |(2 1 / 3 — 
1) = 0.1949 . . .. We remark that Soundararajan's method does not allow us to replace the fi-symbol in 
this assertion by + or fi_. 

1.3. The mean square in the circle problem. The results reported so far already suggest that 

inf{A : P(x) «a x x } = \ . (1.12) 
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To prove (or disprove) this conjecture is usually called the Gaussian circle problem in the strict sense. In 
favor of this hypothesis, there are quite precise mean-square asymptotic formulas of the shape 

x 

J(P(x)) 2 dx = CX 3 / 2 + Q(X), (1.13) 
o 

with 

C = 3^ E r\n)n-^ = £ ^^(1 + « 1.69396 . (1.14) 

n=l 

The estimation of this new remainder term Q(X) has been subject of intensive research, using more 
and more inge nuity. We mention the results of H. Cramer [THj: Q(X) = 0(X 5 / 4+t ), E. Landau [TPS] : 
Q(X) = 0(X 1+£ ), A. Walfisz HZ5|: QpT) = 0(X(logX) 3 ), and I. Katai US]: Q(X) = 0{X{\ogX) 2 ). 
Later on, E. Preissmann |143| found a short and elegant proof for this last result, using a deep inequality 
of H.L. Montgomery and R.C. Vaughan [Tl7| . 

Recently W.G. Nowak |138j succeeded in improving Katai's estimate to 

Q(X) = O (x(logX) 3/2 log 2 x) . (1.15) 

We shall sketch the idea of this refinement. On any interval [hX, X], one can split up P(x) as 

P{x) = H(X,x)+R(X,x) , 

where 

x 

J R(X, xf dx < X 1/2 , 

X/2 

and 

j H(X, X fd X = c(x^-(lX)^) +0 (x £ ^S))' (1-16) 

with 

Ab(«) := min \k — n\ , B :— {n £ Z : r(n) > } . 

This is deduced by standard techniques, applying (1.6) and a Hilbert-type inequality of H.L. Montgomery 
and R.C. Vaughan |TT7| . 

Using the trivial bound Ab(h) > 1 in (1.16), along with the well-known fact that ^2 n ^ y r(n) 2 -C 
ylogy, one readily obtains the result of I. Katai in a way that mimicks E. Preissmann's |143| 

argument. The basic idea of the improvement in (1.15) is the observation that the set B is actually 
rather sparse, i.e., Ae(n) should be "on average" somewhat larger. It was known to E. Landau |1U2) that 

1 ~ /i — ( x -> °°) . 
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for some c > 0. Roughly speaking, this suggests that the elements of B behave like some sequence 
[d n Vlog n] . For the latter, the distance between subsequent terms asymptotically equals d Vl°g n - If 
one could use this bound for Ab(ii) in (1.16), the estimate (1.15) would be immediate. However, the 
following result, which can be proved, and turns out to suffice for this purpose: For every integer h ^ 0, 

£ r{n)2<< (im^ x{logx)1 ' 2 ' (L17) 



\<n<x, n+h£~B 



where <f> denotes Euler's function. This inequality is derived by an elementary convolution argument from 
the related bound 



y i « ( _WM _Z 



1 < 11 < X 



\4>{\h\k) J logx ' 



which is true for arbitrary integers k > and h =/= 0. The latter estimate follows by Selberg's sieve 
method, and was established (for k = 1) by G.J. Rieger |149) . 

We conclude this subsection by the remark that, during the last 15 years, higher power moments of 
P(x) have been investigated as well. We mention the papers of K.M. Tsang |165| . D.R. Heath-Brown 
A. Ivic and P. Sargos [HI], and W. Zhai |182j which provide results of the shape 

(P{x)) m Ax - (~l) m C m X 1+m/4 (X^oo), (1.18) 

where 3 < m < 9 and C m > are explicitly known constants. This not only supports the conjecture 
(1.12), but is of particular interest in the cases of odd m: It means that there is some excess of the 
values x for which P{x) < over those with P(x) > 0. This observation matches well with the different 
accuracy of the and ^--estimates (1.9) and (1.10). 

2. Lattice points in spheres. 

As most experts agree, in our familiar three-dimensional Euclidean space the analogue of the Gaussian 
problem is the most difficult and enigmatic. For integers n > 0, let 

r 3 (n) := #{(m, k, I) € I? : m 2 + k 2 + I 2 = n } . 

Then the formula 

J2 r 3 (n) = ^x 3 / 2 + P 3 (x) (2.1) 

0<n<x ' 

defines the lattice point number, as well as the lattice point discrepancy P 3 (x\ of a compact, origin- 
centered ball of radius y/x. By the same argument which furnished (1.2), it is trivial to see that P?,{x) <C x. 
However, it was known already to E. Landau |l()3j that 

P 3 (x) = O (x 3 / 4 ) . (2.2) 

For many decades, the subsequent history of improvements of this O-bound was essentially the personal 
property of I.M. Vinogradov, whose series of papers on the topic in 1963 culminated in the result |172| 

P 3 (x) =o(x 2 ' 3 {\ogxf^ . (2.3) 
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Also in this case, the last few years brought significant further progress in this problem, namely the 
estimates of F. Chamizo and H. Iwaniec 

P 3 (x) = O (x 29 / 44+e ) (2.4) 

and of D.R. Heath-Brown 

P 3 (x) = O (x 21 / 32+e ) . (2.5) 

Note that || = 0.65909 . . ., || = 0.65625. These bounds are based on the finer algebro-arithmetic theory 
of 7*3(71), notably on the explicit formula, an account on which can be found, e.g., in P.T. Bateman 0. 

In contrast to this quite long history of upper bounds, the sharpest known O-result was obtained by 
G. Szego as early as 1926 |160| and reads 

P 3 (x) = n_ (V^logiz*) 1 / 2 ) . (2.6) 

Neither Hafner's nor Soundararajan's ideas were able to improve upon this bound, probably because the 
function 7*3(71) is distributed quite evenly and is not supported on a small subset of the integers. 

However, the corresponding f2 + -estimates have a longer history of their own: After K. Chandrasekha- 
ran's and R. Narasimhan's result m] that 

I* F 3(X) 

hm sup — = +°° ' 

x — >oo X * 

W.G. Nowak [HHj proved that 

P 3 (x) = n + (x^ 2 (\og 2 xfPQoga x)- 1 ' 2 ) . (2.7) 

Later on, S.D. Adhikari and Y.-F.S. Petermann pP obtained the improvement 

P 3 {x) = n + (x 1/2 \og 2 x) . (2.8) 

Finally, and again quite recently, K.M. Tsang |168| complemented Szego's bound (2.6), showing that in 
fact 

P 3 {x)=Sl ± (x 1 / 2 {\ogx) 1 ' 2 ) . (2.9) 

As for the mean-square of the lattice point discrepancy P 3 {x), its asymptotic behavior has been 
investigated by V. Jarnik as early as in 1940 |S7]. He showed that, with a certain constant C3 > 0, 

x 

J{P 3 {x)) 2 &x = c z X 2 \ogX + o(x 2 {\ogX) 1 l 2 ^ . (2.10) 


The proof of (2.10) is much more difficult than in the planar case. It uses the theory of thcta-functions, 
complex integration and the classic Hardy-Littlewood method. As far as the authors were able to ascer- 
tain, an improvement of the error term in (2.10) has never been attained. 
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The interesting point of (2.10) is that it contains the omega estimates (2.6), resp., (2.9), in full 
accuracy, apart from the ambiguity of sign. That is, our state of knowledge is quantitatively quite 
different than with the circle problem: While (1.13) together with (1.9) - (1.11) tell us (in the planar 
case) that P(x) <C x 1 / 4 in mean-square, with unbounded sequences of "exceptional" x-values for which 
P(x), resp., —P(x), is definitely much larger, a phenomenon of this kind is not known for the three- 
dimensional sphere. 

For spheres of dimensions s > 4, the situation is much easier and better understood, leaving no room 
for progress. With the elementary identity 

r 4 (n):=#{ft,fc 2 ,t3,t 4 )eZ 4 : k{ + k\ + fcjj + kj = n } = 8 a(n) - 32 a ( J) , 

m | n, m>0 

at one's disposal, it is not difficult to show that (see E. Kratzel [EH], p- 227-231) 

2 

P±{x) := r 4 (n) - y x 2 = 0(a; logx) , P 4 (x) = Q (x log 2 x) . (2.11) 

0<n<x 

These rather simple estimates have been improved slightly to 

Pi{x) = O [x{\ogx) 2 l^ , P A (x) = fl± (x \og 2 x) . (2.12) 

The O-result in (2.12) is due to A. Walfisz |178j . |179j who used some extremely deep and involved 
analysis, while the £!-bounds are due to S.D. Adhikari and Y.-F.S. Petermann pQ. 

For dimensions s > 5, the exact order of the lattice point discrepancy P s (x) of an origin-centered 
s-dimensional ball of radius y/x, has been known for a long time, namely 

P a (x) = O (x 3 ' 2 - 1 ^ , P s (x) = Q (x 8 ' 2 - 1 ^ . (2.13) 

See again E. Kratzel [SH], pp 227 ff. For the finer theory of these higher dimensional lattice point 
discrepancies, the reader is referred to the monograph of A. Walfisz |177| . 

3. Lattice points in convex bodies. 

By a p-dimensional convex domain CD P we mean a compact convex set in W with the origin as an 
interior point. Its distance function F is a homogeneous function of degree 1, such that 

CD p = {te R p : F{t- CD p ) < 1}. 

The number of lattice points b£Z p in the dilated convex domain xCD p , where a; is a large positive 
parameter, is denoted by 3 

A(x; CD p ) = S V : F{n; CD p ) < x}. 

3 Observe the difference in notation compared to sections 1 and 2: If CD P is the p-dimensional unit ball, then A(x; CD P ) = 
A p (x 2 ). 
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Let area(C D2) denote the area of CD 2 and vol{CD p ) the volume of CD p (p > 2). The number of lattice 
points in xCD p can be written as 

A(x;CD 2 ) = area(CD 2 )x 2 + P{x;CD 2 ), 
A{x;CD p ) = vol(CD p )x p + P{x;CD p ). 

Here area(CD 2 )x 2 and vol(CD p )x p denote the main terms and P(x;CD p ),p > 2, the remainders or 
error terms. 

3.1 Lattice points in plane convex domains. An elementary result due to M. V. Jarnik and H. 
Steinhaus |159| . 1947: Let J be a closed, rectifiable Jordan curve with area F and length I > 1, and let 
G be the number of lattice points inside and on the curve. Then \G — F\ < I. 

Conclusion. Let I denote the length of the curve of boundary of the convex domain CD 2 . Then 

\P{x;CD 2 )\ < Ix. 

The result P(x; CD 2 ) <C x is the best of its kind under the above conditions. An example is the square 
with length of side 2 and (0, 0) as centre, dilated by x. 

Boundaries with nonzero curvature. Suppose that the boundary of the plane convex domain CD 2 
is sufficiently smooth with finite nonzero curvature throughout. We say that the boundary curve is of 
class C k (k = 2, 3, . . .) if the radius of curvature is k — 2 times continuously differentiable with respect to 
the direction of the tangent vector. Let r be the absolute value of the radius of curvature in a point of 
the boundary and r max ,r m j n the maximum, minimum of r, respectively. Then we shall assume that 

< r min < r max < 00. 

Van der Corput's estimate. The theorem of J. G. van der Corput |17Uj . 1920: Let the curve of 
boundary of the plane convex domain CD 2 be of class C 2 , then 

P(x;CD 2 ) < xi. 

A more precise formulation due to H. Chaix |10|. 1972: 

\P(x;CD 2 )\ < C(l +r max x)3, C > an absolute constant. 

A formulation with explicit constants is due to E. Kratzel 92 , 2004: 

|P(x;Ci?2)|<48(wa ; )* + f703Vw + ^^ N )Vi+ll for x > 1. 

In particular, if CD 2 = E 2 is an ellipse, defined by the quadratic form 

Q(ti,t 2 ;E 2 ) = at 2 1 + 2bt 1 t 2 +ct 2 2 , ac-b 2 = 1, 

then one has 

\P(x; E 2 )\ < 38xi + (iai + 700^) ^ + 11. 
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E. Kratzel and W. G. Nowak 2004, showed that the factor 38 is the leading term of the estimate 
may be reduced to 4r. In particular, for the circle E 2 = C (a = c = 1, 6 = 0), 



\P(x;C)\ < 



^-X3 + f^/i + 3 for x > 28, 
9xi for x > 1000. 



V. Jarnik JBI], 1925, proved that the estimate of J. G. van der Corput is the best of its kind. In fact, 
there exist plane convex domains CD 2 with P(x; CD 2 ) = ul(x 2 / 3 ). 

J. G. van der Corput |171| . 1923, showed that his estimate can be sharpened for a large class of 
plane convex domains C'D 2 . Under the additional assumption that the boundary curve is of class C°°, 
he obtained 

P(x;CD 2 ) <x' 5 

holds with some 1? < f . O. Trifonov |164| . 1988, showed that the estimate holds with = |j + e,e > 0. 



3 . ^. "^uwi liiiil! x-uvyu, uxxvy,»v,vy. uxxv, xxvyxvux, « 41 

Huxley's estimate. The theorem of M. N. Huxley g2], 1993: Let the curve of boundary of the plane 
convex domain CD 2 be of class C 3 . Then 

P(x;CD 2 ) <x^(logx)^, — = 0,6301.... 

M. N. Huxley 2003, improved this result and obtained 

131 18637 131 

P(x;CD 2 ) < x^(logx)^, — = 0,6298.... 

208 

G. Kuba 1994, proved a two-dimensional asymptotic result for a special ellipse E 2 . Let 

{2 2 
(m,n)G Z 2 :^ + ^<1 
x y 

Then it is shown that, for xy — ► oo, 

A(x, y; E 2 ) = nxy + P(x, y; E 2 ) 

with 

' / i 23 n \ Mi 119 , 315 

(xy) 73 (log xy) us for x <C y 100 (log xy) 100 , 



P(x,y:E 2 )<$: 
In the last case one also has 

P(x,y;E 2 ) = n( _ j 



■j= for x 3> y 100 (logxy) wo . 



Lower bounds. Suppose that the boundary curve of the plane convex domain CD 2 is of class C 2 . Then 
the following lower bounds have been proved: 

The Theorem of V. Jarnik (13 . 1924: 

P(x;CD 2 ) = Q fx*) . 
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S. Krupicka [31], 1957: 

P{x;CD 2 ) = fi± fxi 
W. G. Nowak HUI, 1985: 

P{x;CD 2 ) = (a:* (logs)*) . 

In this connection the estimates of the average order of the remainder are interesting: 
S. Krupicka 1957: 

X 

J \P(Vt; CD 2 )\dt > art. 
o 

W. G. Nowak H25|, 1985: 

y |P(Vt;OD 2 )|di<ar4, y |P 2 (Vi;CA>)|dt <xi 



P. Bleher [7], 1992: 

a; 

J P 2 {t;CD 2 )dt~ Ax 2 , 
o 

A is a positive constant depending on CD 2 . 
W. G. Nowak (EH, 2002: 

x+A 



J P 2 (t;CD 2 )dt~4:AAx, 



x—A 

for any A = A (or) satisfying 



z^oo A(arJ 



Points with curvature zero at the boundary. 

Lame's curves. Let k E N, k > 3. A Lame curve L^ is defined by 

|*i|* + M fc = i. 

The curvature of this curve vanishes in the points (0,±1), (±1,0). The number of lattice points inside 
and on the dilated Lame curves xLj- is defined by 

A(x; L k ) = #{(m, n) e Z 2 : \m\ k + \n\ k < x k }. 

The points with curvature zero give an important contribution to the estimate of the number of lattice 
points. 

A first estimate. The number of lattice points is represented by 



11 



2 r 2 (j) 
*r(f) 



Trivially A(x; L k ) <C a;. The first progress concerning the upper bound was made by D. Cauer ||], 1914, 
who proved that 

A(x;L k ) < x 1 "^, 

which was also obtained by J. G. van der Corput |169) . 1919. The resolution of the problem of the size 
of 

A(x;L k ) = 0, 

was given by B. Randol |144) . 1966, for even k > 2 and by E. Kratzel [75], 1967, for odd k > 3 and by 
the same author [75], 1969, for k = 3. 

The second main term. The phenomenon of a second main term was introduced by E. Kratzel |75|. 
1967, EU, 1969 (see also [EH]): 

A(x; L k ) = x 2 + 4*l$ k (x) + P(x; L k ). 

The function (x) is defined for k > 1, and for ^ > i it is represented by the absolutely convergent 
infinite series 

^ + 1 -^)E(^) 2 4 fe) (27rna;), (3.1) 

' n— 1 

where {2tttix) denotes the generalized Bessel function (see E. Kratzel [53], p- 145). It is known that 

i^f k (x) = 0, ^(V-i) . 

Accordingly, this term in the asymptotic representation of A(x; L k ) is called the second main term. The 
remainder P is estimated by 

P(x; L k ) = o (x — * oo) 

for k > 3. 
Improvements: 

E. Kratzel [HO], 1981 and independently W. G. Nowak [EI], 1982: 

27 27 

P{x;L k ) < x^T, — = 0,6585.... 
W. Muller and W. G. Nowak HD], 1988: 

P(x;i fe )«xl(logx)M, g =0,6578.... 

G. Kuba US], 1993: 

P{x;L k ) «x"(logx)^, ^ = 0,6301.... 
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Lower bounds: 

L. Schnabel [HjQ, 1982: 

P(x;L k )=Sl(xty . 

E. Kratzel US], 1988: 

P{x;L k ) = Q ± 

W. G. Nowak H2H, 1997: 

P(a;; J L fe ) = r>_ (^(logx)^ 
M. Kiihleitner, W. G. Nowak, J. Schoissengeier and T. Wooley [27, 1998: 

P(x;L 3 ) = n + (at a (fog log a;)* J . 

The average order: W. G. Nowak [T52|.1996: 



M. Kiihleitner 2000: 



X 

~[ P 2 {t;L k )dt < x. 
o 

x 

~ [ P 2 (t;L k )dt = C k x + {x 1 - s ) 



with explicitly given Cfe>0,e = e(fc)>0. 
W. G. Nowak [Eg, 2000: 

P 2 (i;L fe )di <x(loga;) 2 . 

x-l/2 

M. Kiihleitner and W. G. Nowak |100j . 2001: Let A(x) be an increasing function with 

. / , x , , logx 
Mx) < - and lim -f- = 0. 

2 !E— >OQ A(XJ 

Then 



J P 2 (t;L k )dt~4C k Ax 



x—A 

with the same constant C k as above. 



General boundary curves. Suppose that the boundary curve of the plane convex domain CD2 contains 
a finite number of points with curvature zero. Y. Colin de Verdiere |16j . 1977, showed that 

P(x;CD 2 ) « a; 1 "* 
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holds if k — 2 > 1 is the maximal order of a zero of the curvature on the boundary of CD 2 . Moreover, this 
estimate cannot be improved if the slope of the boundary curve is rational in at least one point, where 
the curvature vanishes to the order k — 2. Further, B. Randol |146| . 1969, Y. Colin de Verdiere |TCj, 1977, 
and M. Tarnopolska- Weiss, 1978, showed: If CD 2 is rotated about the origin by an angle ip then 

for this new domain CD' 2 the estimate 

P(x;CD' 2 ) < x% 

holds for almost all tp € [0, 2ir}. 

W. G. Nowak |125| . 1984, obtained a refinement in the rational case. Assuming that the boundary 
curve has rational slope in an isolated point with curvature zero, then the contribution of this point to 
the discrepancy in the asymptotic representation of the number of lattice points can be given explicitly 
by a Fourier series which is absolutely convergent. E. Kratzel [85], 2000, gives an integral representation 
for this contribution. It has the precise order x 1 ~ 1 ' k if k — 2 is the order of vanishing of the curvature 
in the boundary point. In other words: Each isolated boundary point with curvature zero and rational 
slope produces a new main term in the asymptotic representation of the number of lattice points. 

Further, W. G. Nowak |127| . 1985, obtained a refinement also in the irrational case also. Assume 
that the slope in a boundary point of vanishing curvature is irrational. Under certain assumptions on the 
approximability of this irrational number by rationals, this point gives a contribution to the error term 
of order x® with ■& < |. Improvements of this result were given by W. Miiller and W. G. Nowak, |119| . 
1985. Moreover, M. Peter, |139| . 2000, proved that there is an error term with § < | if and only if the 
slope is of finite type at the point with curvature zero. 

3.2 Lattice points in convex bodies of higher dimensions. An elementary result is the Theorem 
of J. M. Wills, |180| . 1973: Let K be a. p-dimensional, strictly convex body (p > 3) and r the radius of 
its greatest sphere in the interior. Let G be the number of lattice points inside and on the body. Then 

voK K) [l-^ P < G <vol{K) (l + ^) P , 

provided that r > y/p /2. 
Conclusion: 

A{x\ CD P ) = vol(CD p )x p + 0{x p - r ). 

Boundaries with nonzero Gaussian curvature. Suppose that the boundary of the convex domain 
CD p (p > 3) is a smooth (p — l)-dimensional surface with finite nonzero Gaussian curvature throughout. 
Assume that the canonical map, which sends every point of the unit sphere in R p to that point of the 
boundary of the convex domain where the outward normal has the same direction, is one-one and real- 
analytic. The first non-trivial estimates for the remainder P(x; CD p ) were given by E. Landau, 1912 
and 1924, in case of an ellipsoid and by E. Hlawka, 1950, for a general convex body. We begin with the 
general case, and we shall give a short survey of the history of the ellipsoid problem. 

The general case. 

The Theorem of E. Hlawka, jSZj, EH, 1950: 
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Improvements: 



E. Kratzel and W. G. Nowak, [53], 1991: 



P(x;CD p ) < 
E. Kratzel and W. G. Nowak [5J], 1992: 



x p 2+7^+4 for 3 < p < 7 ; 
:r P- 2 +3FFT for 8 < p. 



J X P 2+ 5 P + 2(l 0ga ;)5p + 2 fo r 3<p<6, 

I ar 2 p (log x) p for 7 < p. 

W. Miiller [122, 2000: 

{x 1+ t§+ £ for p = 3, 

x 2 +A+ £ for p = 4, 

p-2+ ,/+ 4 +e r _ . 

X F fo r 5 < p 

for every e > 0. 

W. G. Nowak [EH, [BO], [HI] showed over the years 1985 - 1991 



P{x; CD P ) 



(x(logx)3^j for p = 3, 



0± ^x~(loga;) 2 2 pJ for p > 4. 

A mean square estimate: W. Miiller |121| . 1997: If p > 4, then 

y \P{t;CD p )\ 2 dt^x 2p - 3+e 



for every e > 0. 

Bodies of revolution. Let F(t, z) be a distance function depending on two variables t, z. Then 

Rp = {(ii,t 2 ,...,t p _i,2)el ! ' : F^t 2 1 +t 2 2 + --- + tl_ 1 ,z) < l} 

denotes a p-dimensional body of revolution. F. Chamizo [12], 1998, proved for the number of lattice 
points A(x; R p ) the representation 

A(x; R p ) = vol{R p )x p + P(x; Rp), 
where the remainder is estimated by 

x~ for p = 3, 
P(x;Rp)<gi{ x 3 logx for p = 5, 
r p ~ 2 for p > 5. 
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M. Kiihleitner 2000, proved the lower bound 

P(x;R 3 ) = (a; (log x) i (log log x) ^ e~ A V lo s lo s lo s , 
where A is a positive constant. 

M. Kiihleitner and W. G. Nowak [TUT] . 2004, improved this to 

P(x;R 3 ) = fi_ (a;(log x)i (log log x)^^ 2 ' 1 ^ (log log logs) - § 

An effective estimate in case of p = 3: E. Kratzel and W. G. Nowak, 2004, proved a qualitatively 
weaker estimate, but with precise numerical constants. We describe only the simplest case. Let the body 
of revolution be given by 

R 3 = {{h,t 2 ,z) G K 3 : t\ + t\ < f(z), \z\ < 1} , 

where it is assumed that 

f(z)>0, f(z) = f(-z), /(1) = 0, /(0)<1, 

/ (z) — ► — oo, / (z) — ► — co for z — > 1 — 0. 

Further assume that f"(z) is monotonic and < r m j n < r max < co for the absolute value of the radius of 
curvature. Let 

A(x;R 3 ) = #{(n 1 ,n 2 ,7n) G Z 3 : nf + n 2 2 < x 2 f 2 (^) , \m\ < x} 
= vol(R 3 )x 3 + P(x;R 3 ). 

Then 

|P(z;-R 3 )| < (73r r J lax + 13r,i ax J x 2 + cixlogx + c 2 x + c 3 x* +c 4 i J +5 
with explicit positive constants ci, . . . , C4 depending on r m i n and r max . 

Ellipsoids. The square of the distance function of an ellipsoid CD P = E p (p > 3) is given by the positive 
definite quadratic form 

p p 

F 2 (t; E p ) — 22 o-v^vifi., a un — a^v G K, 

V— 1 fl— 1 

with the determinant c? = det{a vp ) > 0. Then the ellipsoid i? p is defined by 

-Ep = {te R p : F(t; E p ) < l} , 

and the number of lattice points is given by 

A(x; E p ) = # {n G Z p : F{n; E p ) < x} 
irP/ 2 



r(| + i)V5 
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x p + P(x;E p ). 



It is necessary to distinguish between rational and irrational ellipsoids. An ellipsoid is called rational 
if there exists a number c > such that ca vil are integers for all v 1 /i. An irrational ellipsoid is then a 
non-rational ellipsoid. 

The Theorem of E. Landau Q3H|, 1915, [W|, 1924: The estimates 

P(x;E v ) = { ^} p-1 , > 



I x~ 

hold for arbitrary ellipsoids. 

Rational ellipsoids. A. Walfisz [T731, 1924 ; for P > 8 and E - Landau [TUSl, 1924 ; for P > 4: 

„ . f x p ~ 2 for p > 4, 

P(x; E p ) < <^ 2 . 2 , ^ , 
v p ' \ x 2 log x for p = 4. 

Note that 

P(x; E p ) = n (x p - 2 ) for p > 3. 

A. Walfisz [T7HI, I960: 

P(x;E 4 ) = O (x 2 (logx)^ 

Y.-K. Lau and K. Tsang [11111, 2002: 



P(x;E 3 ) = n ± (xy/k>gx). 



Irrational ellipsoids. In the case of irrational ellipsoids a lot of special results are available. To obtain 
them, one usually imposes special conditions on the coefficients of the quadratic form. Most of them were 
proved by B. Divis, V. Jarnik, B. Novak and A. Walfisz. It is nearly impossible to list all the different 
results. For details, see F. Fricker's monograph |24], 

The most important results have been obtained for quadratic forms of diagonal type: 

p 

F 2 {t;E p ) = ^a„ v tl cw>0. 



V. Jarnik ^3, 1928: The estimate 

P{x; E p ) « 



x p 2 for p> 4, 
x 2 log 2 x for p = 4 



holds for arbitrary diagonal quadratic forms. 
V. Jarnik 1929: The estimate 

P(x; E p ) = o(x p - 2) for p > 4 
holds for irrational, diagonal, quadratic forms. 
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V. Jarnik [SHI, 1929 for k = 4 and A. Walfisz [T75| . 1927 for p > 4: For each mono-tonic function f{x) > 
with f(x) — ► for x — ► oo there is an irrational, diagonal, quadratic form such that 

P(r- FM - J fi (^VW) for p > 4, 

K ' p) ~\Q.{x 2 f(x)\og\ogx) for p = 4. 

V. Bentkus and F. Gotze g], 1997: The estimate 

P(x; E p ) < x p " 2 for p > 9 

holds for arbitrary ellipsoids. This result was extended to p > 5 by F. Gotze : 26 , 2004. 
V. Bentkus and F. Gotze [5], 1999: The estimate 

P(x;E p ) = o{x p - 2 ) for p>9 

holds for irrational ellipsoids. 

Points with Gaussian curvature zero at the boundary. 

Super Spheres. Let k,p £ N, k > 3, p > 3. The super spheres SSk, p are defined by 

|t 1 | fe + |*2| fe + -- + |^| fe = l. 

The Gaussian curvature of these super spheres vanishes for t v — [y = 1, 2, . . . ,p). Consider the number 
of lattice points inside and on the super spheres xSSk.p' 

A(x; SS k , p ) =#{neZ p : K| fc + \n 2 \ k + ■■■ + \n p \ k < x k } . 

A first estimate: B. Randol |145| . 1966, for even k and E. Kratzel [78], 1973, for odd k: 

■2\ p rf(|) 



A(x; SS k , P ) = V k , p x p + A(x; SS k , P ), V k , p = 



kj r(i + |) 



A(x; SS ktP ) 



0(x p ~ 2+ ^} for fc<p+l, 

0,Q (xb-W-ty for k>p+l. 



Hence the lattice point problem is settled in case when k > p + 1, as in the planar case. 
An asymptotic representation: The Theorem of E. Kratzel, 1988: 

p-i 

A(x; SS k , p ) = V k;P x p + ^ H k, P A x ) + Afc, p (x). 

r=l 

The terms Hk lPl r{x) are defined recursively: 

1 



A kil (x) = -2^(x) = -2 [ x - [x] 
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Hk, P ,r(x) = (JJ (p ~ r)^ fclP _ r y - O -1 " 1 ** -1 ^*) dt 

o 

for r=l,2,...,p — 1. Furthermore 

H k ,p,i(x) = pV k , p -iif)^f k {x), 
where ipi k \x) is defined in (3.1). Estimates: 

f x p_2 for p — r > k, 



r = 2, 3, . . . ,p — 1, and 



ff fc> p, r (aO = n± (a;(f- r )( 1 -i)+ r ^) for r = 2, 3, . . . ,p - 1, 

A fc!j3 (x) = ^± • 

The representation of the number of lattice points suggests that besides of the main term Vk lP x p there 
are p — 1 further main terms. But the known upper bounds for these terms at present allow at most one 
second main term. So far we have 

[ V kp x p + (x p - 2+ -^) forfc<jj + l, 

Afar SSi ) = < ^ ' 

' \ + pVfcj.-iV'i/i^) + ° (^ (p - 2)(1 - for k>p + l. 

Improvements: 

E. Kratzel |H3], 1988, for p = 3 and R. Schmidt-R6h [122, 1989 for p > 3: 



Afe, p (or) « ^ p _ 2+ ^_ 



x p~ 2 + Tifi for 3<p<7, 
£P- 2 +3Fft for 8<p. 



Estimates of the third and fourth main terms: E. Kratzel [S3], 1988, for r = 2 and S. Hoppner and E. 
Kratzel [D], 1993, for r = 3: 

X 

A fe , 2 (^)dt < x 1 "^, 
A fci3 (ti)di < xlogx 
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and, if A fc , r (i) < A* r (x) for 1 < t < x, 

H kiP , r (x) « X (P-^ + ^) (AI^X)) 1 - 2 ^ (log x) 



r-2 



provided that r = 2,3, k > p — r. It is highly probable that the estimate also holds for r > 3. 
Conclusion: G. Kuba's estimate for A^ ^(x) and E. Kratzel's estimate for Afc 3 (a;) lead to 

H k , Pt2 (x) « x^-^-^+^ilogx)^ 1 -^ for k>p-2, 
H k , p , 3 (x) < x (p - 3)(1 "* )+ M for £;>p-3. 

In case when r > 4, p > 5 it is true that 

£ff fc , p , r (o ; )«^- 4 )( 1 -ft)+*. 

Summary: 

{ V k . p xP + O {x Xk ") for fe<p + l, 

^' p) = { V k [ p xP +pV k ^ ( p %(x) + O (x*>*) for k > p + 1 



with 



_ f P -2 + ^ i for 3<p< 7, 
Afc ' P " \p-2 + 1 ^ ri for 8<p, 
12 

$k, P = P-2 + -— for p = 3,4, k<p + 4, 
7p + 4 



(p- 


2)(1- 


165 \ 
146k) 


-1- 73 -1- t 


for 


p = 3,4, 


(p- 


2)(1- 




- - + — 

5 T 5fc 


for 


P>5, 


(p- 


2)(1- 


165 \ 
146fe/ 


-1- 73 -1- t 


for 


P>5, 



53^482^ £>0 



168 

In addition to all that, E. Kratzel [SSj, 1999, proved that the estimate holds with 

5 



#fc,p =p-2 



3p+l 

for even k and p + 1 < k < 2p — 4. 

Super ellipsoids. Let k,p £ N, k,p > 3. Consider the super ellipsoids SE P 

Xl\tl\ k + \2\t2\ k + ■ ■ ■ + \\tp\ k = 1 

with Ai, A2, . . . , X P > and the number of lattice points inside and on the super ellipsoids xSE p : 
A{x-SE p )=#{n£V> ■.\ 1 \n 1 \ k + \ 2 \n 2 \ k + --- + \ p \n p \ k <x k ). 



20 



V. Bentkus and F. Gotze 0, 2001, proved 



A(x; SE P ) = (AiAz , P A yj- k * P + SE P ), 

where 

P(x; SE P ) < x p - k 

for even k and sufficiently large p. A super ellipsoid is called rational if there exists a number c > such 
that all c\j are natural numbers. Otherwise the superellipsoid is called irrational. Then 

P(x;SE p ) =o(x p - k ) 

if and only if SE p is irrational. 

General convex bodies. Our knowledge about the properties of points on the boundary with Gaussian 
curvature zero for the estimate of the number of lattice points is very poor. Something is known if the 
points with Gaussian curvature zero are isolated. 

K. Haberland [2ZJ, 1993: Let the boundary of the convex body contain only finitely many points of 
vanishing curvature such that the tangent plane is rational at these points. Then each such point produces 
a new main term. 

E. Kratzel (HOj, 2000, and (EH|, 2000, simplified the proof of this result for p = 3 and gave integral 
representations for these contributions. Furthermore, if the slope of the tangent planes are irrational the 
error terms will be of smaller order. 

The restriction to isolated zeros of curvature excludes some important bodies such as the super spheres, 
for example. Thus M. Peter, |140j . 2002, extended the considerations to the case of non-isolated zeros of 
curvature. Instead of assuming only finitely many zeros of curvature, he assumed only finitely many flat 
points. 

The reader is also referred to the papers of E. Kratzel, :9l , 2002, and D.A. Popov, j!41| . 2000. 
4. Divisor problems and related arithmetic functions 

The classical (or Dirichlet) divisor problem consists of the estimation of the function 

A(x) := d ( n ) ~ z( lo gz + 27 - 1) - |, (4.1) 

where d(n) — 1 i s the number of divisors of a natural number n, 7 is Euler's constant, and the 

prime denotes that the last summand in (4.1) is halved if x ( > 1) is an integer. 

As the reader will observe in the results to come, there is a far-reaching analogy between this error 
term A(x) and P(x), the lattice point discrepancy of the circle, which we discussed in Section 1. The 
deep reason is the great similarity of the generating functions 

E^T=C 2 (s) and jr r M = 4C m (s)=H( S )L(s) (Re(s) > 1) . 
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Here C®(i) is the Dedekind zeta-function of the Gaussian field, and L(s) is the Dirichlet L-series with the 
non-principal character modulo 4. 

The generalized (Dirichlet) divisor problem (or the Piltz divisor problem, as it is also sometimes 
called) , consists of the estimation of the function 

A fe (z) := dk(n) - xP fe _i(logz) - (4.2) 

where dk(n) is the number of ways n may be written as a product of k given factors, so that di(n) = 1 
and d2(n) = d(n), A2(x) = A(x). In (4.2), Pk~i{t) is a suitable polynomial in t of degree fc, and one has 
in fact 

P k -x{\ogx) = Res x s -\ k (s) S -\ (4.3) 

S— 1 

where = Xm>i n ~ S (^ es > 1) * s the Riemann zeta-function. The connection between dk(n) and 
C k (s) is a natural one, since one has 

oo 

C k (s) = Y,Mn)n- s = TJ (l-p-T k (Res>l). (4.4) 

n— 1 p prime 

From (4.4) one infers easily that dk{n) is a multiplicative function of n and that, for a prime p and a a 
natural number, dk(p a ) = k(k + 1) ■ • ■ (k + a — l)/al. The relation (4.4) may be extended to complex k; 
for this and other properties of the so-called generalized divisor problem see |60| . Chapter 14. The basic 
quantities related to Ak(x) are the numbers 



Qk 



:= inf { a k > : A k (x) < x a " }, j3 k := inf { b k > : J A 2 k (x)dx < X bk }. (4.5) 



Starting from the classical result of Dirichlet that A (a;) <^ ^/x, there have been numerous results on ctk 
and Pk', for some of them the reader is referred to [HO], Chapter 13, |162j . Chapter 12, These 
results, in principle, have been obtained by two types of methods. For k = 2, 3 the estimation of ak is 
carried out by means of exponential sums, and for larger k by employing results connected with power 
moments of C( s )- The results on j3k have been obtained by several techniques, also using power moments 
of C(s). When k = 2 an important analytic tool for dealing with A(ir) is the formula of G.F. Vorono'i 
(see e.g., Q]J], [57j, [§U|, [HE!, Q75]). It says that 

OO 

A(x) = -2 7 r- 1 V^5I d ( n ) n ~ 1/2 (Ki(^Vnx) + |}i(47rVnx)) , (4.6) 

n=l 

where Ki,Yi are Bessel functions in standard notation (see e.g., (lllj ). Despite the beauty and importance 
of (4.6), in practice it is usually expedient to replace it by a truncated version, obtained by complex 
integration techniques and the asympotic formulas for the Bessel functions. This is 

A(x) = (ttv^)- V/ 4 ^ d(n)n- 3 / 4 cos(4tt % A^ - tt/4) + O e (x e {\ + {x/N) 1 ' 2 )} , (4.7) 

n<N 

where the implied O-constant depends only on e, and the parameter N satisfies 1 <C N < x A for any 
fixed A > 0. Estimates of the form 

A(x) < x a (logx) c (4.8) 
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with various values of a and C > have appeared over the last hundred years or so and, in general, 
reflect the progress of analytic number theory. The last in a long line of records (see jSU], Chapter 13 
for a discussion) for bounds of the type (4.8) is a = if| = 0.3149 . . . , due to M.N. Huxley This 
result is obtained by the intricate use of exponential sum techniques (see his monograph 02|) connected 
to the Bombieri-Iwaniec method of the estimation of exponential sums (see the works of E. Bombieri 
H. Iwaniec jS] and of H. Iwaniec - C.J. Mozzochi [H2|)- The limit of these methods appears to be 
oil < 5/16, whilst traditionally one conjectures that a 2 = 1/4 holds. In general, one conjectures that 
ctk = (k — l)/(2k) and f3 k — (k — l)/(2k) holds for every k > 2. Either of these statements (see j^O] 
or |162| ) is equivalent to the Lindelof hypothesis that C(| + it) <C e \t\ £ for any given e > 0; note that 
trivially (3 k < ct k holds for k > 2. 

The best known upper bound a.3 < 43/96 was obtained by G. Kolesnik [72], who used a truncated 
formula for A3 (a;), analogous to (4.7), coupled with the estimation of relevant three-dimensional expo- 
nential sums. Further sharpest known bounds for ak, obtained by using power moments for £(s) are (see 
pT| . Chapter 13): a k < (3k - 4)/(4fc) (4 < k < 8), a 9 < 35/54 and (see |E] for the values of k between 
10 and 20), a k < (63/c - 258)/(64fc) (79 < k < 119) (see For large values of k one has the best 

bounds which come from the Vinogradov-Korobov zero-free region for £(s) and the bound for £(1 + it). 
Namely if D > is such a constant for which one has 

aa + it)^\t\ D ^ 3 ' 2 \og 2 l s \t\ (i<(7<l), (4.9) 

then 

a k < l-i-2 2/3 (L>fc)- 2/3 . (4.10) 

From the work of H.-E. Richert |148j it is known that D < 100 may be taken in (4.9). The best known 
result D < 4.45 is due to K. Ford |25]. 

As for results on f] k , it is worth noting that the exact value (3 k = (k — l)/(2fc) is known for k = 2, 3 
(see HOI, Chapter 13) and k = 4 (see E2). One has 5 < 9/20 (see W. Zhang [THU), fa < 1/2, fa < 
0.55469, P 8 < 0.606167, (3 9 < 0.63809, /3 10 < 0.66717 (see \g2\). If (4.9) holds, then (ibid.) 

0k < l-l{Dk)- 2 '\ (4.11) 

On the other hand, one may ask for lower bounds or omega results for A k (x). G.H. Hardy showed in 
[5T] that 

Further improvements are due to K.S. Gangadharan K. Corradi- I. Katai J7| and to J.L. Hafner 
|28j . K. Corradi- I. Katai proved 

A(x) = (x 1/4 exp{c(log 2 x) 1/4 (log 3 a;)" 3 / 4 }) (c> 0) . 

Hafner op. cit. showed that, with suitable constants A > 0, B > 0, 

A(aO = fl+ ((slogx) 1 / 4 (log 2 s)( 3 + lo s 4 )/ 4 exp(-A v /Ioi 3 ^ : )) . (4.12) 
Recently K. Soundararajan |157j (see his lemma in Section 1.2) proved 

A(x) ^^((xlogx) 1 / 4 ^^)' 3 / 4 ^ 2473 - 1 )^^)- 5 / 8 ) . (4.13) 
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Note that (3/4)(2 4 / 3 - 1) = 1.1398 , whilst (3 + log4)/4 = 1.0695 , but on the other hand 

(4.13) is an omega result and not an Q + or f2_ -result. In other words, Soundararajan obtains a better 
power of log 2 x, but cannot ascertain whether A(x) takes large positive, or large negative values. His 
method, like Hafner's (see HHI), carries over to Ak(x). Hafner proved, with suitable > 0, 

A, (a:) = n* ((xlogx)( fe - 1 )/( 2fe )(log 2a; )« fe - 1 )/ 2fe » felo s fe - fe + 1 )+ fc - 1 exp(-A fc Vlog^)) , 

where f2* = fi + if k = 2,3 and CI* — Sl± if k > 4. Soundararajan showed that 

A k (x) = fi ((xlogx)( fe - 1 )/( 2fe )(log 2 x)« fe+1 )/ 2fe )( fc2fe/<fc+1, - 1 )(log 3 a ; )- 1 / 2 -( fe - 1 )/ 4 ' £ ) . 

The above estimate holds with fi + in place of £1 if k = 3 (mod 8), and with in place of il if k = 
7 (mod 8). 

Large values and power moments of A(x) were investigated by A. I vie It was shown there that 



/ \A(x)\ A dx < e x 1+A ' A+e (0<A<35/4), 
Jo 



(4.14) 



and the range for A for which the above bounds holds can be slightly extended by using newer bounds for 
C,{\ + it). The bounds in (4.14) are precisely what one expects to get if the conjectural bound a.2 = 1/4 
holds. They were used by D.R. Heath-Brown to prove that the function x~ 1 ^ 4 A(x) has a distribution 
function and that, for A £ [0, 9] (not necessarily an integer), the mean value 



X- l - A ' A [ X \A(x)\ A dx 
Jo 



converges to a finite limit as X — > oo. Moreover the same is true for the odd moments 

r-X 



X- 1 -^ 4 f A(x) A dx (A = 1,3,5,7,9). 
Jo 

For particular cases sharper results on moments are, in fact, known. A classical result of G.F. Voronoi 
|173j states that 

/ A(x)dx = -X + 0{X 3 ' 4 ), 
so that A(x) has 1/4 as mean value. K.C. Tong |H)3) proved the mean square formula 

/ A 2 (a;)dx = CX 3/2 + F(X) (C = (Btt 2 )" 1 ^ d 2 (n)n~ 3/2 = 0.6542869 .. .), 

Jl n=l 

where the error term F(X) satisfies F(X) <C ^log 5 X. Much later E. Preissmann |143| improved this to 
F(X) <§; Xlog 4 X. Further relevant results on F(X) are to be found in the works of Y.-K. Lau- K.M. 
Tsang |109j and K.-M. Tsang |165) . |166j . In particular, it was shown that 

f x 

J F(x)dx = -i&ir 2 )-^ 2 log 2 X + cX 2 \ogX + O(X) (4.15) 
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holds with a suitable c > 0. On the basis of (4.15) it is plausible to conjecture that, with a certain k, one 
has 

F(x) = -(47r 2 ) _1 a;log 2 x + Kxlogx + 0(x). 
K.-M. Tsang |167| treated the third and fourth moment of A(x), proving 



J A 3 (x)dx = BX 7 ^ + O e {X l3+e ) (B>0) 



(4.16) 



and 

rX 



J A 4 (a;)dx = CX 2 + £ (X 7+E ) (C > 0) (4.17) 

with f3 = §g>7 = §§• In a forthcoming work of A. Ivic and P. Sargos jSJ, those values are improved to 
f3 = |,7 = || . A result on integrals of A 3 (a;) and A 4 (x) in short intervals, which improves a result of 
W.G. Nowak |137| . is also obtained in The analogues of (4.16) and (4.17) hold for the moments of 
P(x), with the same values of the exponents (3 and 7. One of the ingredients used therein for the proof of 
(4.17) is the following lemma, due to O. Robert- P. Sargos |15()| : Let k > 2 be a fixed integer and S > 
be given. Then the number of integers ni, 112,113, 77,4 such that N < nx,ri2, 713,714 < 2N and 

\n\ ,k + n\ /k -n\ /k -n 1 J k \<5N 1 ' k 

is, for any given e > 0, 

< £ N e {N 4 5 + N 2 ). (4.18) 

Moments of A k (x) when 5 < k < 9 are treated by W. Zhai 182], who obtained asymptotic formulas with 
error terms. 

Finally we mention results for the Laplace transform of A(x) and P(x), studied in [23 an d We 
have 

J P 2 (x)e- X / T dx=^(^-) ^r 2 (n)n- 3 / 2 -T + O e (T Q + £ ) (4.19) 

and 

/•OO -1 / rp \ 3/2 OQ 

/ A 2 {x) c - X / T dx = - — ) V rf 2 (n)^ 3 / 2 + r(A 1 iog 2 r + A 2 iogr + A 3 ) + Oe(r' 3+£ ). (4.20) 

The Aj's are suitable constants (Ai = — l/(47r 2 )), and the constants \ < a < 1 and 5 < /3 < 1 arc 
defined by the asymptotic formula 

^ r(n)r(n + /*) = V ^ ^(-l) d d + E(x, h), E(x, h) < £ x a+e , (4.21) 
d{n)d(n + h) = x Y,(^g xY ^ ( l ^) J +D( ^ x ' h ^ D{x, h) ^ £ x^ . (4.22) 

n<a i=0 j=0 d\h ^ ' 

The Cy's are certain absolute constants, and the ^-bounds both in (4.21) and in (4.22) should hold 
uniformly in h for 1 < h < x 1 ! 2 . With the values a = 5/6 of D. Ismoilov 03| and (5 = 2/3 of Y. Motohashi 
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[TT5| it followed then that (4.19) and (4.20) hold with a = 5/6 and (3 = 2/3. Motohashi's fundamental 
paper (op. cit.) used the powerful methods of spectral theory of the non-Euclidean Laplacian. A variant 
of this approach was used recently by T. Meurman |116) to sharpen Motohashi's bound for D(x, h) for 
'large' h, specifically for a; 7 / 6 < h < x 2 ~ e , but the limit of both methods is (3 — 2/3 in (4.22). Using the 
results of F. Chamizo ^3], A. Ivic obtained a = 2/3 in (4.19), which appears to be the limit of the 
present methods. 

The general divisor problem can be defined in various ways. Here we shall follow the notation intro- 
duced in [STj. Let d(a%, a 2 , ■ ■ ■ akl n ) be the number of representations of an integer n > 1 in the form 
ri = m^m^ 2 . . • m^ k , where k > 2 and 1 < <Xi < 02 < ■ • • < a k are given integers, and the m's are positive 
integers. Then by the general divisor problem we shall mean the estimation of the quantity 

k / k \ 

A(ai,a 2 , • ■ -a k ;n) = ^ d(ai, 02, • • • afc! n) - ^ Res I C(a,«) J a; s s 1 

k / k \ 

= ^d(ai,a 2 ,...a fe ;n)-^ [] CKMOU 1 ^ (4.23) 

if the Oj's are distinct. If this is not so, then the appropriate limit has to be taken in the above sum. For 
instance, if a\ = a% = . . . = = 1, then we get the classical Dirichlet divisor problem (without the dash 
in the sum and the constant term in (4.2)). It should be noted that k does not have to be finite. A good 
example for this case is 

00 

C(s)C(2s)C(3s)... = < n ) n ~ S (Res>l), (4.24) 

n=l 

corresponding to k = 00, aj = j for every j. Using the product representation (4.3) (with k = 1) for C( s )j 
it follows that a(n) is a multiplicative function of n, and that for every prime p and every natural number 
a one has a(p a ) — P(a), where P(a) is the number of (unrestricted) partitions of a. Thus a(n) denotes 
the number of nonisomorphic Abelian (commutative) groups with n elements (see |82| for a survey of 
results up to 1982). 

Another example of this kind is 

00 

Y[ C{rm 2 s) = J2s(n)n- S (Res>l), 

r>l,m> 1 7i—l 

where S(n) denotes the number of nonisomorphic semisimple rings with n elements (see J. Knopfmacher 

eh). 

In the general divisor problem we introduce two constants: 

a k :— a(ai, a 2 , ■ ■ ■ , a k ) = inf { a > : A(ai, a%, . . . , o/t) < x a } , 
Pk~P(ai,aa,...,Ok) = in{l/3>0:j A 2 {a u a 2 , . . . , a k )dx « X 1+2(i \ , 
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which generalize a k , (3 k in (4-5). From the classical results of E. Landau it follows that 

k- 1 k- 1 f \ 

<a k < - ■ - (4.25) 



2(oi + a 2 + . . . + a k ) (fe + l)oi 

if the numbers dj are distinct. In many particular cases the bounds in (4.25) may be superseded by the 
use of various exponential sum techniques, coupled with complex integration methods and the results on 
C(s) (moments, functional equations etc.). The case k = 3 (the three-dimensional divisor problem) is 
extensively discussed in E. Kratzel's monograph jHH] and e.g., the work of H.-Q. Liu |113j . 

In what concerns results on A. Ivic [SJ proved the following result: let r be the largest integer 
which satisfies (r — 2)a r < ai + a 2 + . . . + a r _i for 2 < r < k, and let 

r - 1 

g k = g( ai ,a 2 ,. . . ,a k ) := Tf ■ ■ ■ r. (4.26) 

2(ai + a 2 + . . . + a r ) 



Then (3 k > g kl and if 

i-T 



|C(| + ^)| 2fe " 2 dt < e T 1+£ (4.27) 



o 



holds, then (3 k = g k . Thus assuming (4.27) we obtain a precise evaluation of (3 k , but it should be remarked 
that (4.27) at present is known to hold for k = 2,3, whilst its truth for every k is equivalent to the Lindelof 
hypothesis. Moreover, we have 

J A 2 (a 1 ,a 2 ,..., afc )d a; = r!(x 1+29 Mog A x), (4.28) 

where A = A(a\, a 2 , . . . , a k ) is explicitly evaluated in |51j . The results remain valid if k = oo, provided 
that the generating function is of the form 

C 91 (MC 92 (MC 93 (M--- (1 <h <b 2 <b 3 < ...), (4.29) 

where the 6's and the q's are given natural numbers. The generating function (4.24) of a(n) is clearly of 
the form (4.29). In this case it is convenient to define the error term as 

6 / oo \ 

n <x j=i V e=i,e^£j ^ J 

Then (4.28) implies that 

J E 2 (x)dx = n{X 4/3 logX), (4.31) 

and (4.31) yields further E(x) = ^(a; 1 / 6 log 1 ^ 2 x), obtained by R. Balasubramanian-K. Ramacha- 
ndra 2 . The omega result (4.31) is well in tune with the upper bound 

J E 2 {x)dx = 0{X 4 / 3 \og S9 X) 
of D.R. Heath-Brown 33 , who improved the exponent 39/29 = 1.344827 ... of gUJ. 
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As for upper bounds for E(x) of the form E{x) — 0(x c \og C x) (C > 0) or E{x) — O e (x c+£ ), they have 
a long and rich history (note that c > 1/6 must hold by (4.31)). The bound c < 1/2 was obtained by P. 
Erdos-G. Szekeres who were the first to consider the function a(n) and the so-called powerful numbers 
(see @H], 123) f° r an account). After their work, the value of c was decreased many times (chronologically 
in the works (7D|, E7J, USD, Q2E1, EH| and Liu QI2] , who had c < 50/199 = 0.25125...). Recently 

O. Robert-P. Sargos |150| obtained c < 1/4 + e by the use of (4.18). This is the limit of the method, and 
the result comes quite close to the conjecture of H.-E. Richert |147j . made in 1952, that E(x) = o(x 1 ^ 4 ) 
as x — > oo. The most optimistic conjecture is that c — 1/6 + e holds. 

A related problem is the estimation of T(x) = X) T (C)i where t(G) denotes the number of direct 
factors of a finite Abelian group G, and summation is over all Abelian groups whose orders do not exceed 
x. It is known (see fRj^ or |S3]) that 

oo 

T(x) = *(«), E *W n ~ S = C 2 (s)C 2 (2s)C 2 (3s) ... (Res > 1), 

n<x n—1 

which is of the form (4.29) with qj = 2, bj — j. Here one should define appropriately the error term (note 
that t(n) = ^2d\n a (d)a(n/d)) as 



Ai(x) := T{x)-^{D j \ogx + E j )x 1 



li 



3=1 

5 

= E a{m)a{n) -^(Djlogx + E 3 )x 1/j , 

mn<.x j— 1 

where Dj and Ej (D% > 0) can be explicitly evaluated (see H. Menzer-R. Seibold |114| 1. They proved 
that Ai(x) = O e (x p+e ) with p = 45/109 = 0.412844 . . . , improving a result of E. Kratzel HU, who had 
the exponent 5/12. As for the true order of Ai(x), one expects p — 1/4 to hold. This is supported by 
the estimate (see [S3] and [M]) 



J &\{x)&x = Sl{xV 2 \og i X) 



(4.32) 



which yields p > 1/4. The value of p was later improved (see |115| and |113p . and currently the best 
bound p < 47/130 is due to J. Wu |181j . Moreover, in [S3] it was shown that 

J* A 2 1 (x)dx = O e (X s / 5+£ ), 

while if (4.27) holds with k = 5, then one has 

J Al(x)dx = £ {xV 2+£ ). (4.33) 

In this case (4.32) and (4.33) determine fairly closely the true order of the mean square of Ai(x). 

The previous discussion centered on global problems involving arithmetic functions which arise in 
connection with divisor problems, that is, the estimation of the error terms in the asymptotic formulas 
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for summatory functions or the related power moment estimates. One can, of course, treat also local 
problems as well, namely problems involving pointwise estimation, distribution of values, and other 
arithmetic properties. The literature on this subject, which can be also considered as a part of the theory 
of divisor problems, is indeed vast, especially on d(n) and dfe(n). Thus it would be a great increase in 
the length of this work, as well somewhat outside the mainstream of the theory, if we dwelt in general 
on local problems. Also one sometimes, under divisor problems, considers divisor problems of the type 
(4.29) (or (4.23)), where the q's are allowed to be negative. This greatly increases the class of functions 
that are allowed (like e.g., squarefree numbers whose characteristic function is generated by £(s)/£(2s), 
or squarefull numbers whose characteristic function is generated by C(2s)C(3s)/C(6s)). In this work we 
have found it best to adhere to divisor problems of the form (4.23) or (4.29). 

As for local problems, we shall conclude with just a few words on a representative divisor problem, 
namely the function a(n). It was proved by E. Kratzel [77] that 

hmsup log a(n) ■ — = , (4.34) 

n^oo log n 4 

and his result was sharpened by W. Schwarz - E. Wirsing |155) and (on the Riemann Hypothesis) by 
J.-L. Nicolas |123| . General results for multiplicative functions, analogous to (4.34), were obtained by 
E. Heppner and A. A. Drozdova - G.A. Freiman (see |142j . Chapter 12). The existence of the local 
densities of a(n), the numbers 

rffe := lim x^ 1 y 1 

x — >oo ^ — J 

n<_x,a(n)—k 

for any given integer k > 1, was proved by Kendall- Rankin |7()j . This was sharpened in |46| to 

l=d k x + 0(x 1/2 \ogx), 

n<-x,a(n)—k 

where the O-constant is uniform in k. Further sharpenings and generalizations are to be found in 47 , 
[S3, EH, E7J and [T55| . 

The distribution of values of a(n) was investigated in 07] and P. Erdos - A. Ivic [20]. It was shown 
there that [x, 2x] contains at least ^fx integers n for which a(n + 1) = a(n + 2) = . . . = a(n + k) with k = 
[logcclogg cc/(40(log 2 x) 2 )}, and at least y/x integers m such that the values a(m + l),a(m+2), . . . , a(m+t) 
are all distinct, where for a suitable C > 

t = [C(\ogx/ log logx) 1/2 ]. 

Iterates of a(n) were treated by P. Erdos - A. Ivic [21] and 
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